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ON . We study the kinetics of bimolecular, catalytically-activated reactions (CARs) in 

d-dimensions. The elementary reaction act between reactants takes place only when 
these meet in the vicinity of a catalytic site; such sites are assumed to be immobile 
and randomly distributed in space. For CARs we develop a kinetic formalism, based 
on Collins-Kimball-type ideas; within this formalism we obtain explicit expressions 
for the effective reaction rates and for the decay of the reactants' concentrations. 
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O ! I- INTRODUCTION 

o 

f^. ■ Many industrial and technological processes depend on catalytically-activated reactions 

(CAR), whose understanding, however, used to be mainly phenomenological 1 . More micro- 
scopically inclined approaches emerged only during the last decades. Thus much progress 
was made in determining how reactions are promoted by specific catalytic substrates 2 . From 
the point of view of many-body effects much understanding was gained from an extensive 
study of the CO-oxidation on metal surfaces 3 . The first type of research 2 is at the molecular 
level and, ideally, leads to K eh the reaction rate for two molecules which meet at a catalytic 
site. However, the results of Refs.3 demonstrate amply that the mere knowledge of K e i is not 
sufficient for detemining the time evolution of the global reaction process. Thus the decay 
forms for the reactants' concentrations in CARs often deviate strongly from the results of 
formal-kinetics schemes 3 . We note that such deviations are not encountered in CARs only, 
but are widespread 4-6 , being due to many-body effects, often associated with fluctuations 
in the spatial distributions of the reacting species. Consequently, a robust, reliable descrip- 
tion of CARs has to go beyond formal-kinetic considerations and to take the influence of 
higher-order particle correlations into account. Here we develop such an approach. 

We start from the following system: The catalytic substrate consists of N, immobile cat- 
alytic sites (CSs), randomly placed in the <i-dimensional reaction volume V. The reacting 
particles of type A are also randomly positioned; their initial average density is no = ua(0). 
The A-particles diffuse with diffusion coefficient D. For simplicity we neglect sterical hin- 
drances between the As and also between As and CSs. Furthermore, we assume that the A 
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particles do not get trapped by the CS to form metastable A-CS complexes. Next, pairs of A 
particles may react upon encounter, if they meet at a CS (Fig.l). Formally, the elementary 
reaction step is 



where K e i(r), is a random function of the spatial variable r; K e i(r) is strictly zero everywhere, 
except in the vicinity of any of the CSs, in which case K e i(f) = K e \. Here we take the CSs 
to be randomly distributed in V with a constant average density ric'- thus the CSs' set and 
V have the same topological dimension. We hasten to add, however, that our approach can 
be readily extended to other distributions of CSs (i.e. regular or strongly inhomogeneous) . 

Since in V many microscopic reactions which follow Eq.(l) go on at the same time, 
the global picture becomes complex, and is, in general, properly described by an infinite 
hierarchy of coupled, differential equations 3,4 . The analytical approach for CARs which we 
pursue here is based on the truncation of this coupled system at the level of third-order 
correlations; the latter are then accounted for through appropriate boundary conditions. 
For non-catalytic reactions in 3d this method was pioneered by Collins and Kimball 8 (CK); 
it generalizes Smoluchowski's treatment of bimolecular chemical reactions (see, e.g. Ref.6). 
The CK-approach provides both for the binary reaction A+A — > 0, and for the recombination 
reaction A + B — > with n^(0) 7^ n B (0) a reasonable description 8 ' 9 of the experimentally 
observed kinetics over the entire time domain. We note, however, that the reaction A+B — > 
with ua(0) = ub(0) requires to go beyond the CK-approach, in order to depict the kinetic 
behavior at very long times, where many-particle effects come into play 4-6,10 . 

Two remarks are here appropriate: 

(a) In the case when the CSs cover V completely, CARs behave exactly as non-catalytic 
reactions. This special case is equivalent to the original CK-problem. 

(b) Recently, the kinetics of the A + B + C — > diffusion-limited reaction has been 
analysed 11 ' 12 using an extension of the Smoluchowski approach (see, e.g. Ref.6). Setting 
A = B = C leads formally to the reaction A + A + A — > 0, for which the procedure shows 
fair agreement 13 between the analytical predictions and the numerical data. For CARs one 
may now be tempted to follow a similar course, by setting A = B in A + B + C — > 0, and 
identifying the Cs with the CSs. We will show, however, that such a procedure does not 
describe CARs well; in Id and in 2d it does not lead to the proper long-time decay and in 
3d it accounts only partially for the effect of CSs. 

The paper is structured as follows: In Section II we formulate the model, by writing 
its basic equations, which allows us to extend the 3d CK-approach to arbitrary d and to 
CARs. Here we also point out the relation between the CK-CARs kinetics in <i-dimensions 
and the trapping problem in (d + <i)-dimensions. In Section III we present explicit results 
for CARs' effective rates and for the reactants' decay; we compare these to the findings for 
non-catalytic reactions. Finally, we conclude with a summary of results in Section IV. 



We formulate our model on a rf-dimensional (d = 1,2,3) lattice with lattice spacing a. 
To each site of the lattice, whose position is specified by the vector r, we assign a time- 
independent variable ricif) which assumes two possible values, namely or 1, depending 
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on whether the site is catalytic, ncif) = 1, or not, ncif) = 0. The catalytic substrate is 
the set of all N CSs; we denote it by {Rk}, where Rk is the vector of the k-th CS. Here we 
take the Rk to be random, independent, uniformly distributed variables. The CSs density 
is ric = V~ x *52r n c(f)i where the sum runs over all lattice sites. The case when {Rk} forms 
different types of ordered geometrical arrays will be discussed elsewhere. 

We start at t — with randomly distributed, identical A particles, with mean density 
tiq. Each A particle moves by jumps to nearest-neighboring sites, the average time interval 
between successive jumps being r. We disregard any excluded volume interactions; thus 
all A particles perform independent random walks, with the associated diffusion coefficient 
D = a 2 /2dr. 

Now, whenever an A particle lands on a catalytic site which is already occupied by 
another A, the two As may react at a rate K e \. Reacting As are immediately removed from 
the system, whereas the corresponding CS remains unaffected. On the other hand, As never 
react at non-catalytic sites. 



A. Evolution of the local density and of the two-point joint density functions. 

Let n(f, t) denote the local density of the As. In continuous-time n(f, t) obeys: 

(2) 



^n(f,t) = A? n(f, *) - K el n c (r) n(r 1: r 2 ;t) 



r\=ri=T 



Here n(fi, r 2 ; t) is the two-point joint density function (i.e. the probability of having at time 
t an A at r-y and another A at f 2 ), and the symbol A^ stands for the following difference 
operator acting on r: 

A P n(f,t) = -2dn(f,t) + $^ra(f',t) (3) 

r',nn 

Here the sum runs over nearest neighbors only, i.e. |f ' — r\ = a. 

Apart from the factor nc{r) stemming from the CS, Eq.(2) is the conventional rate 
equation for A + A — > (see e.g. Ref.8 for a discussion). On the rhs of Eq.(2) the first 
term accounts for the particles' migration, while the second one describes the reaction: In 
standard fashion, the reaction term is taken to equal the product of the rate K e i, of the 
probability of having a pair of A particles at the same site and at the same time, and of 
ncif); the latter factor is new here and is due to the fact that As can react only at catalytic 
sites. We note that Eq.(2) embodies mean-field assumptions: In a more rigorous approach 
one has to use the three-body probability that two As encounter each other at a CS; this 
probability is here decoupled by having it represented as the product of ricif) and n(f, r; t). 

Before we turn to the analysis of the time evolution of n(fi,f2',t), it is instructive to 
consider the result of the simple kinetic approach, which later will serve as reference. Taking 
the volume average of both sides of Eq. (2) and assuming that the non-linear reaction term 
on the rhs of Eq.(2) decomposes into a product of averaged local densities, we arrive at the 
standard, formal-kinetic "law of mass action" (see, e.g. Ref.7) 

j t n A {t) = -K el n c n 2 A {t) (4) 
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Here riA(t) = V 1 J2r n (?, t) is the average density of the As. From Eq.(4) the decay of riA(t) 
at sufficiently long times, t ^> l/ncK e in(0), follows: 

n A (t) « (n c K el t)- x (5) 

It is important to note that n A {t) given by Eq.(5) is independent of the spatial dimension 
and is the same in, say, Id and 3d. Second, the formal-kinetic approach predicts that the 
effective reaction rate constant is ncK e i, and is thus independent of other parameters, such 
as, e.g. the particles' diffusivities. 

We now turn to the analysis of the time evolution of n(fi,f 2 ;£). We find that this 
quantity obeys (see, e.g. Ref.13) 

^n(fi,f 2 ;t) = { A n + A ft } n(fi,f 2 ;f) + T, (6) 

in which the terms in curly brackets stem from the particles' motion, while T is a combination 
of joint three-point distributions and arises due to the reaction between the As. 

Eqs.(2) and (6) are the first two equations of an infinite hierarchy of coupled differential 
equations (CDE). Such a hierarchy of CDE cannot be solved exactly; in order to compute the 
evolution of n(f,t) and thus of its mean value riA(t) one has to resort to some approximate 
methods. 



B. The Collins-Kimball's 3d problem and its extension to arbitrary d. 

Here we continue by recalling Collins and Kimball's (CK) analysis 8 of reaction kinetics: 
As discussed, their problem is identical to 3d CARs in a "completely" catalytic medium, i.e. 
such that n c = 1. In Ref.8 the hierarchy of CDE is truncated at the level of the third-order 
correlations, i.e. T in Eq.(6) is set to zero, and the reaction between particles is accounted 
for by introducing a "radiation" boundary condition on n(fi,f 2 ;£); CK stipulate that the 
local reaction rate at any point r, i.e. K e \ n(f,f;t), should be exactly equal to the diffusive 
current of pairs of A particles into this point. Accompanied by the following boundary and 
initial conditions: 

™(ri,f 2 ; Hoo -> n 2 A (t), (7) 

and 

n(fi,f 2 ;t = 0) = n 2 A (0) = njj, (8) 

which signify that correlations in the particles' positions vanish at large separations and that 
initially the particles are uniformly distributed in V, the CK-approach results in a closed 
system of linear equations; these allow then to calculate the effective reaction rate and the 
time evolution of n^t). We note that in their original paper 8 CK have only considered 
the 3d case. Clearly, however, in the formulation of Eq.(6), extending the CK-approach to 
systems of arbitrary spatial dimension d is straightforward, as we show in the following. 

Solving Eq.(6) (with T = 0) subject to the radiation boundary condition and the condi- 
tions in Eqs.(7) and (8), CK find 8 that n A (t) obeys in 3d: 
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j t n A (t) = - K eff (t) n 2 A (t) (9) 

Eq.(9) is similar to the "law of mass action" displayed in Eq.(4), since for both the rhs 
are proportional to the second power of n A (t), in agreement with the elementary reaction 
act. The difference between Eqs.(9) and (4) is that now K e ff(t) stands for an effective (in 
general, time-dependent) rate coefficient, which arises mainly from the particles' diffusion. 
Now K e ff(t) can be expressed through its Laplace-transform: 

K eff {X) = / dt exp(-At) K eff {t), (10) 

J 

for which CK find, for identical reacting particles in the limit A <C D/tq, where r denotes 
the reaction radius: 

A + A, (id 



K eff (X) K el 8nDr 

The second term on the rhs of Eq.(ll), i.e. X/8nDr Q , stems from the Laplace-transformed 
Smoluchowski constant; the latter equals the particles' current towards the surface of a 
single, immobile, adsorbing 3-dimensional sphere of radius r . 

This result is actually a special case of the general, d- dimensional expression for K e ff(X), 
which can be derived within the framework of the CK approach, as we briefly outline now. 

Consider the behavior of the two-point joint density function n(fi, r 2 ; t). Within the CK 
approach it obeys Eq.(6) with T = 0: 

|U(ri,r 2 ;t) = D {A* + A f J n(ri, f 2 ; t), (12) 

where A^ denote the continuous-space Laplace operators. Note that Eq.(12) is a (d + d)- 
dimensional diffusion equation. Now, Eq.(12) is to be solved subject to the initial and 
boundary conditions of Eqs.(7) and (8), as well as to the CK radiation boundary condition 

K el n(r u r 2 -t) = D {V* + V r -J n(n, f 2 ; t)^^ (13) 

As already observed by CK for non-catalytic reactions, n(fi,f^;t) depends only on z = 
\ri — r 2 |, the relative distance between r\ and r 2 , and consequently, the (d + d)-dimensional 
Eqs.(12) and (13) reduce to the effectively <i-dimensional equations: 

Q d 2 d - 1 d 

- nW) = 2D{- 2+ —-}n(z-,t), (14) 



and 

dn(z; t) 



K el n(z; t) = 2 D a d 



dz 



(15) 



z=r 



Here o d = 2ir d ^ 2 rQ~ 1 T~ 1 (d /2) denotes the surface area of the <i-dimensional sphere of radius 
r , with T(x) being the Gamma-function |fL4| . Equations (14) and (15) generalize to arbitrary 
dimensions the 3d equations studied by CK 8 . 
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Now, the solution of Eq.(14) subject to the conditions of Eqs.(7), (8) and (15) can be 
readily obtained by Laplace transforming with respect to t. The Green's function G^(z; A) 
of the equation 



A G d (z; A) = 2 D { + d -± ^} G d (z; A) (16) 



obeys 



A 



= (wW )M/4 ^ ( V^'- (17) 

in which K v {x) is the modified Bessel function 14 . The solution of the Laplace transformed 
Eq.(14), (note the occurence of an additional term, — Uq, on the lhs, as compared to Eq.(16)), 
has the form: 

71 ^ 

n(z;X) = -f [1 + A(X)G d (z;X)}, (18) 

in which the conditions Eqs.(7) and (8) have already been incorporated. In Eq.(18) A(X) is 
a constant, which has to be chosen in such a way that the boundary condition, Eq.(15), is 
also satisfied. Substituting Eq.(18) into Eq.(15) we find that A(X) is: 

A(X) = - G- d \r ;X) {1 - ^ [lnG d (r ; X)]T\ (19) 

-ft el 

where the prime stands for the derivative with respect to r . 

This provides a complete solution of the CK-problem in d- dimensions. In order to cal- 
culate the effective rate constant, K e ff(X), we have merely to insert Eq.(18) with (19) into 
the rhs of Eq.(15). This yields 

Keff(X) = - ^ [lnG d (r )}> {1 - ^ [InG^)]'}- 1 (20) 
Equation (20) takes a physically more revealing form, if we rewrite it as 

1 A + («) 



K eff (X) K el K d (D;Xy 

where 

K d (D; A) = - ^ [InG^ = ^ J^^L (22) 



Eq.(21) resembles the electrostatic law of addition of parallel resistances and displays the 
combined effect of two controlling factors - the effect of the elementary reaction act and the 
effect of the transport of particles towards each other. Notice now that for K d = oo the rate 
K d (D; A) is the Laplace transform of the diffusive current towards the surface of a single 
immobile, perfectly adsorbing <i-dimensional sphere of radius r . Consequently, K d (D; A) is 
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the (i-dimensional analog of the Smoluchowski constant and Eq.(21) represents the desired 
generalization of Eq.(ll) to the d- dimensional case. We hasten to remark that Eq.(21) has 
already been obtained in Ref.10 in a different framework, based on the analysis of the third- 
order joint density functions; Eq.(21) also follows from the general results of Ref. 15, which 
considered reversible reactions. 

The behavior of Kd(D; A) depends in a fundamental way on the spatial dimension d. As 
is well-known (see e.g. Ref.9), in low dimensions (d < 2) Kd(D; A) tends to zero for A — > 
(t — > oo). To be explicit, in the limit A <C D/r^ the parameter K d (D; A) determined by 
Eq.(22) equals (8/J/A) 1 / 2 for d = 1 and 8nD/(X ln(8D/r$\)) for d = 2. In higher dimensions 
(d > 2) K d (D;X) approaches constant values at long times; the values for K d (D;X), for, 
say d — 3, 4 and 6 turn out to be 8nDr /X, 8'k 2 Dtq/X and 871 3 Dtq/X, respectively. The 
dependence of K e ff(X) on d, arising due to the second term in Eq.(21), leads, especially for 
d < 2, to deviations of the decay laws from Eq.(5). 



C. The CK-approach extended to CARs. 

After this overview of the CK approach we now extend it to CARs. Following Ref.8 we 
truncate the hierarchy at the level of the third-order joint density functions, which yields 
the continuous-space Eq.(12). Next, in our system only the encounters which happen on a 
CS may lead to a reaction. Thus, the radiation boundary condition is to be imposed on the 
CSs only. Denoting by Sk the surface of the k-th CS, we hence have that instead of Eq.(13) 
the following boundary condition should be imposed on n(ri,r 2 ;t): 

K e m(f u r 2 ;t)\ r ^ 2€Sk = D {V* + VfJ^.f^)^^ (23) 

For simplicity we take in the continuum Sk to be the surface of the ci-dimensional sphere of 
radius r centered at R^- Equations (12) and (23) are complemented by the initial condition, 
Eq.(8), and by the boundary condition, Eq.(7); this constitutes a closed system of linear 
equations which allows the computation of n(fi, f 2 ; t). In turn, knowing n(fi, r 2 ; t) and using 
Eqs.(2) and (23), we have 



jn A {t) = - ^^{V rl + Vft}n(fi,f 2 ;f) 



(24) 

fi,r 2 eS k 



which defines the evolution of the property of interest, namely of n^(t), the mean density 
of As. 



D. Relation between CARs and the trapping reaction. 

We note now that it is expedient to view Eqs.(12), (23), (7) and (8) from a somewhat 
different perspective, which will allow us to find eventually an approximate analytical solu- 
tion for them, and to explain, on simple physical grounds, some seemingly surprising results. 
As a matter of fact, what the CK approach enables us to do it is to reduce the problem 
of computing the rates of binary reactions, taking place in <i-dimensional catalytic systems, 
to the analysis of (imperfect) trapping in (d + <i)-dimensional systems. The latter problem, 
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and especially its quantum mechanical counterpart, the scattering of quantum particles by 
immobile impurities, have been extensively investigated (see, e.g. Refs. 15-18 and references 
therein) . 

The relation between the CAR and the trapping problem can be most simply illustrated 
for the binary CAR in d — 1. Note that Eqs.(12) and (23) describe the evolution of the 
local density of some compound particles, moving with diffusion coefficient 19 Dona two- 
dimensional plane (ri, r 2 ), where r\ and r 2 are scalar variables. The particles may disappear 
at the locations (R k ,R k ) of the traps, placed on the diagonal r\ = r 2 (Fig.2). Physically, 
each compound particle is a pair of A particles, whose coordinates on the one-dimensional 
line are r\ and r 2 respectively; consequently, r\ = r 2 = R k are the only points where two 
As may enter into reaction, in which case the compound particle may be destroyed by the 
trap at Rk with a finite probability related to K e \. In this language, the reaction rate of the 
CAR, i.e. the rhs of Eq.(24), attains a quite lucid meaning: It equals the volume-averaged 
diffusive current of compound particles through the S k . One can now note the fundamental 
distinction between systems with ric = 1 and with nc < 1: The former case corresponds to 
a situation in which the traps cover the diagonal r\ = r 2 completely; here n(ri,r2) depends 
only on the relative distance between r\ and r 2 , i.e. one has n(ri,r 2 ) = n(r\ — r 2 ), which 
then reduces the problem to Id. On the other hand, in the case tiq < 1 the traps fill only 
some portion of the diagonal, and thus the kinetics remains essentially two-dimensional, 
since n(r 1; r 2 ) depends on both spatial variables; compound particles can cross the diagonal 
harmlessly, i.e. with zero reaction probability, through the gaps between the traps. 

Returning now to the general problem of CARs in <i-dimensional media, the correspond- 
ing mapping leads to considering a trapping problem involving compound particles diffusing 
in a (d + <i)-dimensional space in the presence of imperfect traps placed on a (/-dimensional 
substrate. 

Several analytical approaches have been developed to describe the kinetics of trapping 
in systems with non- uniform spatial distributions of traps (see, e.g. Refs. 6 and 15). We will 
search for the solution of Eqs.(12) and (23) in the spirit of the Green's function method 17 . 
Here we merely outline the steps involved, and address the reader for more details to Refs. 15 
and 17. 

One starts with the Laplace-transformed Green's function solution Gd + d(p; A) of the 
(d + c?)-dimensional diffusion, Eq.(12), 

G -^'^'il» M " 2 ^'/l l * (25) 

in which p = (r*i,r 2 ). Furthermore, the Laplace-transform of n{fi,r 2 ;t) is represented as a 
series in which each term is the Green's function solution of Eq.(12), centered around the 
position of the /c-th trap, 

71 ^ — * 

nirur^X) = -± + £ C k (X) G d+d {p ~ A), (26) 
A k 

with the d + rf-dimensional vector O fc = (R k ,R k ). Eq.(26) obeys automatically Eqs.(12),(7) 
and (8); the coefficients C k (X) are to be choosen in such a way that Eq.(23) is satisfied. 
Substituting Eq.(26) into Eq.(23) we arrive at the following system of iV linear equations 
for the C fc (A): 
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- -o = {i^ 1 + G d+(i (r ; A)} C,-(A) + £' G d+d (®k - A), (27) 

A k 

where j = 1, N, and the prime indicates that the sum in Eq.(27) runs over all k with the 
exception oik — j. 

The exact solution of Eqs.(27) for a given distribution of {Rk} requires the inversion of the 
random matrix \\G d+d (Ok ~ ©j! A)||, see Ref.17. Neglecting fluctuations in the distribution 
of the CSs, in which case the Eqs.(27) simplify considerably, one obtains 

n 2 

C(X) * " \{K^ + G d+d lr ;\) + M scr y (2§) 

where M scr denotes the screening integral (or "shielding" integral in the formulation of 
Ref.15) 



M scr 



< J2 ' Gd+Sk - A) > « n c j j dndf 2 5(n - f 2 ) G d+d (p- A) (29) 



In Eq.(29) the brackets denote averaging over the distribution of {Rk}, and the integrations 
with respect to the variables ri and r 2 extend over the whole volume occupied by CSs, 
excluding the volume of a rf-dimensional sphere of radius r . In the following we turn to the 
limit N, V — > oo, while keeping the ratio N/V fixed, iV/V = n C - We note that Eq.(29) is only 
approximate, since excluded-volume aspects between the CSs are neglected; this limits the 
applicability of the expression to tiq sufficiently small. We note also that within our CK-type 
description the dependence of the effective reaction rate on the geometry of the catalytic 
substrate and/or the distribution of the CSs enters only through the screening integral 
M scr . Consequently, any other geometry of the substrate (it can be, for instance, a two- 
dimensional convoluted surface of porous materials, imperfect crystallites with broken faces, 
kinks and steps, or polymers in solution 1 ) can be accounted for by the use of the appropriate 
distribution functions and by corresponding integrations in Eqs.(29). In particular, the 
details of the averaging procedure in the case when the integrations extend over Gaussian 
polymer chains in solution have been discussed in Ref.20. 

III. RESULTS. 

Now, combining Eqs.(26),(28) and (29), we find from Eq.(24) that n^t) obeys the effec- 
tive "law of mass action" in Eq.(9). In the limit A A^ = D/r$, (t ^> td = r^/D), 
when M scr and G d+d (r ; A) reduce to M scr w n c (A K d (D; A)) -1 and G d+d (r ; A) w 
(A K d+d (D/2; A)) -1 , the effective rate constant attains for ucTq 1 the form 

1 A 1 1 

+ y/n .u + Z r, TTtToTTT ' (3°) 



K eff (X) n c K el K d {D-\) n c K d+d (D /2;A)' 

which represents the desired generalization of the CK-type result, Eq.(21), to CARs. Equa- 
tion (30) is the main result of our analysis and allows to compute n^t), which is related to 
K e ff(t) through Eq.(9). Hence: 
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n A {t) = n {l + n f dr ^//(r)}" 1 (31) 

«/ 

Let us consider first the 3d case. We recall the explicit forms of the parameters Kd(D; A), 
presented in the text after Eq.(22), so that in 3d Eq.(30) takes the following form (nc?"o < 1): 

A A A (32) 



K eff (X) n c K el 8nDr A^Dr^nc 

which signifies that in the limit t — > oo the effective rate constant K e ff(t) approaches a 
constant value: 

1 +^ + ^" (33) 



K ef f n c K e i 8nDr 4n 3 Dr$n c 

Consequently, from Eq.(31) we have for riAit) in three-dimensions and large t: 

n A (t) » (K eff t)-\ (34) 

where K e ff is given by Eq.(33). Equation (34) signifies that in 3d catalytic systems riA(t) is 
inversely proportional to t, i.e. is qualitatively the same as the formal-kinetic Eq.(5). Now, 
Eq.(33) differs from Eq.(5), since here K e ff replaces ncK d . Note from Eq.(33) that K e ff 
reduces to ncK e i only for D — > oo. Thus for finite D, the rate K e ff depends both on the 
reaction radius and on the mean density of CSs. The difference between the CK-result for 
^c r o ~ 1 an d f° r CARs with ucTq <C 1 is the last term in Eq.(30), which is due to the 
Green's function solution of the diffusion equation in 6d. For small values of the parameter 
ucTq and for Dr\ <C K e i, the last term in Eq.(30) provides the dominant contribution to the 
effective rate constant K e ff. In this case Eq.(30) reduces to the result of Ref.12, obtained 
for the three-body problem. 

Consider next the evolution of n^(t) in low dimensional systems, i.e. for d — 1 and d = 2. 
We have from Eq.(30) that in Id K e ff(X) is given by (nc^o < 1): 

A , ( A 1/2 , \ln{AD/\rl) 



K eff (X) n c K el K 8D J AttDuc { J 

Now, Eq.(35) shows that in Id catalytic systems the kinetics is richer than in the 3d case: 
comparing the different terms in Eq.(35) one readily notices that depending on A each of 
these terms may dominate K e ff(X); hence a succession of different kinetic regimes may be 
observed in the time domain. When 4nD/K e i is sufficiently large and nc-r is sufficiently 
small, so that n c r <^ exp(—4irD/K e i), we can consider three different intervals, namely 

X D exp(-47iD/K el ) < A < X D , (36) 

Ad (n c r ) 2 < A < X D exp(-AixD / K el ), (37) 

and 

A < A D (n c r ) 2 (38) 
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(a) In the regime described by Eq.(36) the main contribution to K e ff(X) comes from the 
first term on the rhs of Eq.(35); in this kinetically-controlled regime K e ff(t) = ncK e i and 
hence coincides with the result of the formal-kinetic approach. This behavior persists until 
tkc ~ td exp(AirD/K e i), which can be rather large for D 3> K el . On the other hand, such a 
behavior is unobservable for K e \ 4nD, exemplified by instantaneous reactions in Eq.(l). 

(b) In the regime described by Eq.(37) K e ff(X) is dominated by the third term on the 
rhs of Eq.(35), so that K e ff(t) ~ 47iDnc/ln(At/TD). This Id CAR expression is reminiscent 
(apart of the factor ric and the replacement D — > D/2) of the classical result of Ref.9 for 
the kinetics of binary reactions in 2d. One thus expects that at this stage n^(t) ~ ln(t)/t. 
The appearence of such an effectively 2d regime for CAR in Id constitutes the principal 
difference between the CAR kinetics for nc^o ^ 1 and the kinetics of noncatalytic binary 
reactions in Id. We also note that such an effectively 2d behavior was predicted in Ref.12 
for Id reactions ofA + B + C^O type as the final kinetic stage; in our case, when C 
does not disappear in the reaction act, Eq.(37) determines only an intermediate transient 
stage, which may be observed for times smaller than a typical time t tr , where t tr = 1/Dn 2 c . 
Consequently, for Id systems with very low densities of catalytic sites (nc^o <^ 1) such an 
effectively 2d behavior can last over extended time periods. 

(c) Finally, in the limit of very small A, Eq.(38), K e ff(X) is determined by the second 
term on the rhs of Eq.(35), i.e. K eff (X) w (8D/A) 1 / 2 . Thus for large t the rate K eff (t) « 
(^D/nt) 1 / 2 is independent of n C - Hence for large times the qualitative decay behavior is the 
same for CARs and for non-catalytic reactions. Actually, this seemingly surprising behavior 
has already been observed numerically 21 and has a simple physical interpretation, which we 
will discuss below. Explicitly, we find that the mean density of A particles decays in the 
limit t ^> t tr as 

n A {t) ~ (n/32Dt) 1/2 (39) 
We turn next to the analysis of CAR in 2d. Now Eq.(30) reads: 

1 _ X + XJn(8D m+ A 2 (4Q) 



K eff (X) n c K el 8nD 4ir 2 Dr 2 n c 

On comparing different terms on the rhs of Eq.(40) we infer that in 2d two different kinetic 
stages may take place. Namely, for A from the interval 



X D exp{ — ^ « X « A D , (41) 

n c K d irn c r$ 



the sum of the first and of the third term, which both have the same A-dependence, deter- 
mines K e ff(X), while for smaller A, such that 



. 8nD 1 . 

A « X D exp( g , 42 

n c K e i im c n 



the main contribution to the effective rate constant is given by the second term on the rhs 
of Eq.(40). 

(a) In the regime described by Eq.(41) K e ff(X) is controlled by the constraints imposed 
by the elementary reaction act and by the diffusion in Ad. For such a regime we obtain 
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nM " ( h + (nctrl (43) 

According to Eq.(41), the regime described by Eq.(43) is a transient one and persists until 
ttr « t d exp((D/n c K el ) + (l/rgra c )). 

(b) In the final stage described by Eq.(42), the main contribution to K e ff(X) comes from 
the second term in Eq.(40). Thus, similarly to the behavior in Id, K e ff(X) in the limit A — > 
is independent of n c . This implies that also in 2d for large times the kinetics for CARs is 
the same as for non-catalytic reactions. We find here 

MO ~ *g^. (44, 

which coincides with the result of Ref.9 obtained for the long-time kinetics of non-catalytic 
reactions. We furthermore note that the very long time decay behavior is reached much 
more slowly in 2d than in Id] the crossover time t tr is in 2d an exponential function of n^ 1 
and is thus substantially larger than its Id counterpart, which goes as n^ 2 . 

Lastly, we discuss the physical origin of the fact that in low dimensions the CARs 
long-time decay is independent of n C - Here, the analogy between the CAR kinetics in 
d-dimensions and the trapping problem in (d + <i)-dimensions again turns out to be very 
fruitful. Let us consider first the case d — 1. The equivalent problem is a d = 2 system 
with compound particles diffusing in the presence of traps placed on the diagonal r± = r 2 , 
(see Fig. 2). Now, it is well-known 9 that even in the presence of a single trap the parti- 
cles' density profile around the trap is not stationary; there is a zone around the trap, 
which is depleted of particles and whose size grows with time as \fi. For the situation 
depicted in Fig. 2 at short times (such that ricVDt <C 1) the depletion zones of different 
traps are well separated from each other: Consequently, at short times the traps act inde- 
pendently and the Laplace-transformed particle current towards a given trap is given by 
Eq.(21), J(A) « K el K 2 (D/2; \)/{K el + K 2 (D/2; A)). The effective rate, which is given by 
Eqs.(23) and (24), is now K e ff(X) = n c J(\); the CAR decay shows in this time-domain an 
effectively two-dimensional behavior going as n^t) ~ ncln(t)/t. This behavior also shows 
up in the results of Ref.12. At longer times, the depletion zones of different traps start to 
overlap and the probability n(r\, r 2 ) of finding a compound particle on the diagonal r\ = r 2 
decreases substantially even in the gaps between the traps. This results in a situation in 
which the array of traps distributed on the diagonal r\ = r 2 acts as an adsorbing line. It 
is not surprising then that the current of particles per trap attains a one-dimensional form 
J(A) ~ K d=1 (D; \)/n c , and that even for n c r C 1 at very long times the evolution of 
riAit) proceeds essentially in the same fashion as for non-catalytic reactions. Similarly, for 
2d CARs we have to analyse the kinetics of trapping in Ad, the traps being now distributed 
on a 2d plane. As before, one can now distinguish between two different temporal regimes: 
one finds first a situation in which the traps act independently, which gives rise to an effec- 
tively Ad behavior 12 . This regime crosses over to a stage at long times, where the depletion 
zones of the different traps overlap, so that the array of traps acts as an effectively 2d ad- 
sorbing plane; the decay of n^(i) is then given by Eq.(41), and is independent of ric- We 
note finally, that such a peculiarity of the trapping kinetics in low dimensional system, asso- 
ciated with the formation of non-stationary depletion zones around traps, has already been 
demonstrated in Ref.15. In particular, it was shown 15 that for 2d systems in which N traps 
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are located inside a circular area of radius R one finds two different temporal regimes: In the 
first, intermediate time regime all traps act independently and the effective rate constant is 
proportional to N. This regime crosses over into a long-time stage, in which the depletion 
zones around different traps overlap; at this stage an array of traps acts as a single trap of 
radius R and the effective rate constant is independent of N. 

IV. CONCLUSIONS. 

We now conclude with the discussion of the obtained results. We find for 3d CARs that 
riA{t) decreases inversely proportional to t, Eq.(34), which agrees with the formal- kinetic 
picture, Eq.(5). Distinct from it, the effective rate constant is less than the formal- kinetic 
value, ricK e i, and depends both on the particles' diffusion constant D and on the reaction 
radius r Q . We also note that in the diluted case, when n^r^ <C 1, the result in Eq.(33) 
shows the same dependence on the system's parameters as the one predicted in Ref.12. The 
very long-time behavior in low-dimensions is somewhat surprising, since the decay turns out 
to be essentially independent of ric- The approach to this asymptotic domain is, however, 
very slow for low densities of CSs, n c rQ <C 1, and thus different decay forms appear at 
intermediate times. The crossover times t tr may be very large, since we find that t tr ~ 1/n^ 
in Id and ln(t tr ) ~ i n 2d For n^r^ <C 1 and for extended period of time the decay 

laws obey in Id and 2d riA(t) ~ ncln(t)/t and ~ nc/t respectively. 
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Figure Captions 



Fig.l. Reactions in catalytic media. Open circles denote immobile catalytic sites; the 
filled circles stand for diffusive A particles. Case (1) shows a situation in which an encounter 
of A particles does not lead to reaction, while in the case (2) the reaction may take place. 

Fig. 2. Open circles on the r\ and r 2 -axis denote the catalytic sites; filled circles give the 
corresponding positions of traps and small black circles denote diffusive compound particles. 



15 





o . o 






• 


• 


(2) 






€) 




o 




(1) 




O 


• o o 






• 


• 


O 


o • 


o 



Fig. 1. Oshanin et al., JCP 




Fig.2. Oshanin et al., JCP 



